These MATLAB scripts were used to generate the coordinates of DC axons and DR axons in an FEM model of SCS (see Howell et al., 2014). The following is a brief summary of the scripts:
NoR, node of Ranvier
MYSA, proximal paranodal segment
FLUT, distal paranodal segment
STIN, internodal segments

The following run for a fiber diameter of 9 microns (um):
(use mrg_arclengths to generate files for other diameters)

dorsalcolumn_coord, coordiantes of dorsal column (DC) axons
- requires two files: patient_cordgeom.txt, axonarc_[fiber diameter]um.txt

dorsalroot_coord, coordinates of dorsal root (DR) axons
- requires same two files as above

mrgaxon_arclengths, arc lengths of first 8 elements of model axon
- NoR-MYSA-FLUT-STIN-STIN-STIN-FLUT-MYSA-...repeats to last NoR
- (ses McIntyre et al., 2002)
- the relationship between internodal length (INL), the length between two successive NoR, and fiber diameter, D, is 100D in this study. This is a reasonable approximation; however, I recommend using the values in the McIntyre et al. work to to determine the nonlinear relationship between INL and D (i.e., INL=f(D)). An 6th or 7th order polynomial works well.

seed_dclayer, called by dorsalcolumn_coord and dorsalroot_coord
- calculates xy positions of axons within dorsomedial aspect of spinal cord, known as the DC area

plot_dcol_droot, used to visualize populations in spinal cord



